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Abstract 

LiNio.5Mn1.5O4 , which has a spinel framework structure, is a promising candidate for the cathode 
material of next-generation lithium-ion batteries with high energy density. We investigate the 
structural transition in LijjNio.5Mn1.5O4 (0 < x < 1) through first-principles calculations using the 
projector augmented wave method with the generalized gradient approximation. We calculate all 
the unique Li-site occupation configurations in a unit cell to obtain the total energies and the most 
stable structures for various compositions. Thermodynamic analysis shows that Lio.5Nio.5Mn1.5O4 
with x = 0.5 is the only stable phase for < x < 1. The decomposition energy is lower than 0.1 eV 
for < x < 0.5, but is distinctly higher for 0.5 < x < 1. The decomposition energy reaches 0.39 
eV at x = 0.75. The ratios of the structures at room temperature are calculated from Boltzmann 
factors by using the energy differences between structures. The crystal structure of the unit cell 
changes gradually from x = to 0.5, but changes markedly from x = 0.5 to 1. This first-principles 
study provides a general evaluation of the variation in the crystal structure with the composition 
of the bulk material, which affects the cyclability of the electrode. 

PACS numbers: 61.50.Ah, 64.60.-i, 64.70.qd, 71.15.Nc 

Keywords: ab initio calculation, thermodynamic stability, crystal structure, electrochemical potential, 



1 



I. INTRODUCTION 



the high energy density of lithium-ion batteries makes them particularly suitable energy 
storage devices for the many new portable electronic devices that are coming onto the mar- 
ket. LiCo02 is widely used as a cathode material for lithium- ion batteries, because it has a 
voltage of up to 4 V versus metallic lithium 1 . Voltage is a fundamental property of batteries; 
higher voltages are closely related to larger electrical energy capacities and faster electronic 
performance in processing units. Recently, several transition-metal-doped spinel cathode 
materials with voltages of 5 V, such as LiM 3 .Mn 2 _ ;r 04 (M = Ni, Cr, Co, Fe, Cu), have been 
reported 2-5 . LiNio.5Mn1.5O4 is a particularly promising material because of its high voltage, 
electrochemical stability, and lack of trace metal impurities 6 ' 7 . Furthermore, materials such 
as Li[MnNiCo]0 2 5 ' 8 , Li[MnNiCr]0 2 9 , and Li[MnNiFeCo]0 2 10 are based on LiNi . 5 Mn 1 . 5 O4 
. The physical and chemical properties of LiNio.sMn! 5 4 have been experimentally de- 
termined by X-ray diffraction, scanning electron microscopy, infrared- Raman spectroscopy, 
electron diffraction, and cyclic voltammetry 11-13 . Although the crystallographic and elec- 
trochemical properties of LiNio.sMni 5 4 have been well reported, the microscopic changes 
in the crystal structure and electrochemical properties during the lithiation and delithia- 
tion of an LL r Nio.5Mni. 5 04 crystal are not yet fully understood. Charging and discharging 
are primary functions of secondary batteries; therefore, the change in the composition and 
structure of electrodes is unavoidable, and the changes influence cyclability. Predicting these 
structural changes will help to accelerate the development of electrode materials. 

Theoretical and microstructural analysis can provide a valuable alternative to experimen- 
tal studies for understanding the properties of materials. First-principles calculation based 
on density functional theory (DFT) is among the most powerful methods for determining 
the fundamental characteristics of electrodes. The theoretical band structures and crystal 
structures of various Li-transition metal oxides, including LiMn0 2 14 ' 15 and LiMn 2 4 16 , have 
been thoroughly investigated. Several studies using theoretical calculations and basic ther- 
modynamic models have shown that the average electrochemical potential of LiNio.5Mnx.5O4 
is derived from the difference in total energy between the fully lithiated and delithiated 
materials 14,17,18 . However, a method has not yet been established for evaluating the changes 
in crystal structure and electrochemical potential with varying Li composition. To analyze 
the properties of crystals at intermediate compositions quantitatively, we investigated the 
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sublattice structures of Li and used DFT to evaluate their thermodynamic stability at var- 
ious compositions between Nio.5Mn1.5O4 and LiNio.5Mn1.5O4 . We propose a method for 
predicting the average electrochemical potential as well as the change in electrochemical 
potential as the composition of the bulk material is varied. This information is useful for de- 
signing new electrode materials, because the electrochemical potential and cyclability of an 
electrode material is strongly affected by the change in the atomic structure during lithiation 
and delithiation. 

II. COMPUTATIONAL PROCEDURE 
A. Crystal Structure 

The starting structure for our study was a unit cell of LiNio.5Mn1.5O4 crystal (Fig. 1). 
The cell belongs to the P4 3 32 space group and consists of 8 formula units containing 56 
atoms in total: LigNi4Mni2032 11 . The O atoms are located at the 8(c) and 24(e) sites, 
and form a distorted face-centered cubic (fee) structure in the cell. Cations occupy the 
ordered polyhedral sites formed by the O sublattice; the Li, Ni, and Mn atoms are present 
in the octahedral 8(c) sites, 4(b) sites, and 12(d) sites, respectively. The Li sublattice forms 
a slightly distorted diamond structure. The eight Li sites in the unit cell were labeled as 
shown in Table I. The initial cells for geometry optimization at various compositions were 
made by removing n (n — 0, 1, 2, • • • , 8) Li atoms from the Li§Ni4Mni2032 cell, resulting 
in a Li 8 _ n Ni 4 Mni2 032 cell. The Li sites are symmetrically equivalent in the space group. 
Therefore, for n — 0, 1, 7, and 8, the unit cell crystal structure is the same. For 2 < n < 6, 
the configurations of the occupied and unoccupied Li sites were considered. We calculated 
all the unique configurations in the symmetry operations of its space group (Table II). The 
geometry of each cell was optimized to obtain the most stable lattice volume, shape, and 
atomic coordinates; the total energies of the cells were computed at the same time. We 
examined differences in energy of the various structures. The change in Gibbs energy is 
expressed by AG = AE + PAV - TAS. However, AE was of the order of 0.04-4 (eV/cell), 
PAV was of the order of 10~ 5 (eV/cell), and TAS* was of the order of the thermal energy in 
this study 14 ' 19,20 . Because the PAV and TAS terms have little effect on AG, the assumption 
that AG can be approximated by only the change in the internal energy, AE, is valid. 
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FIG. 1. (Color online) Crystal structure of LiNio.5Mn1.5O4 . The Li, Ni, Mn, and O ions are 
shown in green, gray, purple, and red, respectively. The radii are different to allow them to be 
clearly distinguished: ru < < r^ n < ro- (a) O polyhedra including cations. The tetrahedra 
(green), light gray octahedra, and dark gray octahedra (purple) indicate the Li04 , Ni06 , and 
MnOg polyhedra, respectively, (b) Li-Li bonding in the Li sublattice. 

TABLE I. Labels and internal coordinates of Li sites in the unit cell. These coordinates are adjusted 
to those of the ideal diamond structure to allow for easy recognition. 

axis L000 L220 Llll L331 L202 L022 L313 L133 
a 0.5 0.25 0.75 0.5 0.75 0.25 
b 0.5 0.25 0.75 0.5 0.25 0.75 
c 0.25 0.25 0.5 0.5 0.75 0.75 



B. Thermodynamic analysis of phases 

The ratio of structures at equilibrium was determined for a total Li content in the bulk 
of x = b. The lithiated and delithiated material was assumed to be Li x X , where X was 
Nio.5Mn1.5O4 . Because this is a pseudo-unitary system of Li, we examine Li content only. 
For compositions of x = a, b, and c (a < b < c) in Li x X , the decomposition reaction is 

U b X '—^U a X + ^^Li c X (1) 
c — a c — a 



TABLE II. Unique configurations of occupied and unoccupied Li sites in a unit celi of 
Li„Ni4Mni2032 . All configurations are assigned an uppercase character from A to W. The filled 
circles and dashes indicate occupied and unoccupied sites, respectively. The labels at the top, for 
example L000, are site names, which are the same as those in Table I. 
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The decomposition energy of this reaction, AGb(a,c), is 

AG b (a,c) = G(b) - (—G(a) + —G(c)) , 
\c — a c — a J 

where G(x) is the Gibbs energy of the most stable structure at composition x. Positive values 
of AGb(a, c) indicate that the left side of Eq. (1) is thermodynamically unstable, whereas 
negative values indicate that it is stable. AG fe (a, c) gives the ratio of the equilibrium amounts 
of species on each side of Eq. (1) by using the Boltzmann factor, 

Pb(a, c) ( AG b (a,c) \ 

where P&(6) and Pb(a, c) are the ratios of both edges in Eq. (1) in the bulk, k is the Boltzmann 
constant, and T is the temperature. T was assumed to be room temperature, 298 K. The 
sum of the ratios of all states at a given a given value of b must be 1: 

P b (b)+ E E (P b (a,c)) = l. (3) 

0<a<6 6<c<l 

The composition was treated as a discrete value. Each ratio of the decomposed state in the 
bulk, Pb(a,c) and P b {b), was calculated by solving simultaneous equations 2 and 3. Then, 
each ratio, Pb(x), was derived from 

&<c<iSn(i,c) (x<b) 

Pb(x) = I P b (b) (x = b) (4) 

^o<a< b ^Pb(a,x) (x>b). 

Assuming that ideal monovalent Li + cations are the only charge carriers, the average 
electrochemical potential versus Li/Li + between x\ < x < x 2 of Li x X , V(x), is given by 

where F is the Faraday constant, and Gu is the chemical potential of Li in a metallic Li 

anode 14 ' 19 ' 21 . 



C. Computational conditions 

The spin-polarized calculations were performed by DFT 22 ' 23 using the plane-wave pro- 
jector augmented-wave method 24 as implemented in the VASP code. 25-27 The exchange- 
correlation term was treated with the Perdew-Burke-Ernzerhof functional 28 . For the unit 
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cell of LizNio.5Mn1.5O4 , a cutoff energy of 500 eV and a /c-mesh of 4 x 4 x 4 were determined 
through preliminary test calculations. A cutoff energy of 500 eV and a /c-mesh of 64 x 64 x 64 
were calculated for Li metal in a primitive cell with body-centered cubic structure. Geome- 
try optimization was truncated when the residual forces on the atoms became less than 0.02 
eV/A. 

III. RESULTS AND DISCUSSION 

A. Thermodynamic stability and phase transition 

Fig. 2 shows the calculated energies of all the structures after geometry optimization. 
For clarity, the values of energies shown on the vertical axis, G'(x), were corrected with 
linearly approximated chemical potential of Li, /i' u : 



The line connecting the energies of the most stable structures for neighboring compositions 
was concave except for the area around x = 0.5 (Fig. 2). This indicates that only one com- 
position, x = 0.5, was thermodynamically stable, except for the end member compositions, 
x = and 1. The Li sublattice in the cell labeled O forms an fcc-like structure. The fee 
structure has high symmetry and uniform distance between atoms. The electrostatic repul- 
sion between Li + ions makes it preferable for them to be positioned as far apart as possible. 
Thus, the high stability of the cell is reasonable. In the region < x < 0.5, the energy 
profile was concave, suggesting that Li x Ni .5Mn 15 O4 should decompose to Nio.5Mnx.5O4 and 
Lio.5Nio.5Mnx.5O4 at K. However, the decomposition energy, AE^O, 0.5), was less than 
0.1 eV. This slight instability would produce the structures labeled B, C, and H at interme- 
diate compositions through the thermal effect, even at room temperature. The most stable 
structure changes were A, B, C, H, and O for compositions varying from x = to x = 0.5. 
Table II shows that the change in the most stable structure between neighboring composi- 
tions can be represented by replacing occupied Li sites with unoccupied sites, or unoccupied 
sites with occupied sites. For example, the only difference between structures H and O is 
at the L022 site. This simple mechanism would result in smooth changes in the atomic 
structure as the composition changes and may contribute to good cyclability. In contrast, 



// Li = G{1) - G(0) 
G'(x) = G(x)-(G(0) + xfi' Li ). 



(6) 
(7) 
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FIG. 2. Energies of crystal structures with a composition of LizNio.5Mn1.5O4 . Each energy on 
the vertical axis is corrected by linear approximation of the Li chemical potential (Eqs. 6 and 7). 
Labels A to W correspond to those in Table II. The solid line connects the most stable structures 
at the various compositions. The dotted line connects all the thermodynamically stable structures 
between Nio.5Mn1.5O4 and LiNio.5Mn1.5O4 . 

the decomposition energy, AE X (0.5, 1), was high in the region 0.5 < x < 1. In particular, 
AE 75 (0.5, 1) was 0.39 eV for structure S. The bulk material should be a mixture of two 
separate phases, namely, structures O and W. Therefore, the crystal structure of the unit 
cell changed considerably in this region. The cells that transformed from Lio.5Nio.5Mn1.5O4 
to LiNio.5Mn1.5O4 were stressed by the volume change, which was estimated to be about 
4%. 

The most stable structures, which are on the solid line in Fig. 2, are dominant at certain 
concentrations because of the Boltzmann factor. Therefore, we can assume that they are the 
only structures in the bulk material. We used the energies in Fig. 2 and Eq. 4 to calculate 
the ratio of structures in bulk LL r Ni .5Mn 1 . 5 O4 at 298 K (Fig. 3). This figure shows that 
the unit cell structure changed gradually via intermediate concentrations when the total Li 
content x was varied from to 0.5, whereas the structure changed markedly from x = 0.5 
to 1. These results suggest that it would be difficult to detect a difference in the internal 
structure of LijNio.5Mn1.5O4 (0 < x < 0.5) by experimental methods. The cyclability of an 
electrode is expected to depend on how smoothly its crystal structure changes. Therefore, it 
is valuable to describe the change in structure quantitatively as in Fig. 3 in order to predict 
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FIG. 3. Calculated ratio of structures in bulk LixNio.5Mn1.5O4 at 298 K. The value of x on the 
horizontal axis indicates the total Li content in the bulk. Darker areas indicate structures with 
higher Li concentration. The uppercase characters in the legend indicate structures corresponding 
to those in Table II and Fig. 2. The value of x for each of the uppercase characters indicates the 
Li occupancy of the structures. 

the durability of new electrode materials. 



B. Electrochemical potential 

There were three thermodynamically stable compositions, where x = 0,0.5, and 1. The 
electrochemical potential in the region < x < 0.5 and 0.5 < x < 1 was calculated by 
using Eq. (5). Fig. 4 shows the variation of the electrochemical potential with the total Li 
content in a bulk electrode. The calculated values of 3.99 and 3.56 V were lower than the 
experimental values of 4.739 and 4.718 V 11 . However, first-principles calculations using the 
generalized gradient approximation typically give underestimates 14 , and the calculated and 
experimental results showed good agreement. The total change in the theoretical electro- 
chemical potentials for compositions of < x < 1 was calculated as 0.43 V. The changes in 
the electrochemical potential shown in Fig. 4 correspond to the changes in the slope of the 
dotted line in Fig. 2 because the electrochemical potential is reflected in /jl u of the material. 
The small change in electrochemical potential resulted from Li .5Ni .5Mn 1 . 5 O4 being only 
slightly stable, that is, having slightly negative decomposition energy, the small stability 
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FIG. 4. Calculated changes in the electrochemical potential of LizNio.5Mn1.5O4 versus Li/Li + . 

of Lio.5Nio.5Mn1.5O4 . If there were no stable phase in intermediate composition of Li x X , 
the material would theoretically have constant electrochemical potential. However, during 
lithiation and delithiation, an intermediate composition that is highly unstable results in a 
sharp change in atomic structure, which is an obstacle to good cyclability. 



IV. CONCLUSION 

In summary, first-principles calculations revealed the microscopic structure and stability 
of LiNio.sMni 5O4 , which is among the most promising materials for lithium-ion battery 
cathodes. We quantitatively evaluated the energies of all the configurations of Li-site oc- 
cupation for intermediate compositions of LizNio.5Mn1.5O4 (0 < x < 1). In this way, we 
assessed the thermodynamic stability of each structure and the change in the electrochem- 
ical potential with varying total Li content in the bulk. Lio.5Nio.5Mn1 5 4 (x = 0.5) was 
the only stable phase for < x < 1. The decomposition energy was lower than 0.1 eV for 
< x < 0.5, but was high for 0.5 < x < 1. The thermodynamic analysis showed that the 
change in structure with varying Li content was gradual for low Li content compositions, and 
rapid for high Li content compositions. The theoretical change in electrochemical potential 
at x = 0.5 was found to be 0.43 V, which reflects Lio.5Nio.5Mn1.5O4 having lower stability 
than Nio.5Mni.5O4 and LiNio.5Mni.5O4 . We have demonstrated that theoretical calcula- 
tions can reveal precise changes in structure, phase, and electrochemical potential and can 
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explain the cyclability of electrode materials. In particular, visualizing the change in the 
crystal structure as shown in Fig. 3 will be particularly valuable for discovering durable 
electrode materials. 
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